Resistive transition and upper critical field in underdoped YBa 2 Cu 3 06+ x single 

crystals. 
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A superconducting transition in the temperature dependence of the afc-plane resistivity of under- 
doped YBa2Cu3C>6+:E crystals in the range T c < 30 K has been investigated. Unlike the case of 
samples with the optimal level of doping, the transition width increased insignificantly with mag- 
netic field, and in the range T c <13 K it decreased with increasing magnetic field. The transition 
point T C (B) was determined by analyzing the fluctuation conductivity. The curves of Bca(T) mea- 
sured in the region T/T c > 0.1 did not show a tendency to saturation and had a positive second 
derivative everywhere, including the immediate neighborhood of T c . The only difference among the 
curves of B C 2(T) for different crystal states is the scales of T and B, so they can be described in terms 
of a universal function, which fairly closely follows Alexandrov's model of boson superconductivity. 



I. INTRODUCTION 

The nature of high-temperature superconductivity is 
presently one of the most interesting subjects of the solid 
state physics. An important topic of research in this 
field is the temperature dependence of the upper critical 
field £> c2 . In conventional (low-temperature) supercon- 
ductors, in accordance with the BCS theory, the curve 
5 c2 (T) is described by a universal function bgcs{t) m 
terms of reduced variables: the temperature is scaled by 
the zero-field transition temperature, t — T/T c , and the 
magnetic field is scaled by the product of T c and the 
derivative of B c2 {T) at T c : b = B /[T c (-dB c2 / 'dT) T=Tc \. 
U The function bscsit) is linear in the neighborhood of 
T c and saturates to b ~ 0.7 at t — 0. In high-temperature 
superconductors (HTSC) the behavior of B C 2(T) is rad- 
ically different. In Tl 2 Ba 2 Cu0 6 § and Bi 2 Sr 2 CuO y §] 
films, and in Ko.4Bao.6Bi03 single crystals, a posi- 
tive second derivative and a sharp increase in _? c2 (T) at 
low temperature have been detected. Similar properties 
of function _> c2 (T) have been observed in other HTSC 
systems, namely, in YBa 2 (Cui_j,Zn a )306+a; with a criti- 
cal temperature lowered by the strong scattering || and 
Sm1.85Ceo.15CuO.4_j, with n-type conductivity. |7[ 

HTSC is not the only class of materials where the up- 
per critical field does not follow the BCS universal func- 
tion bBcs{t)- But, as concerns HTSC, such deviations 
are probably present in all materials of the family, and 
magnitudes of these deviations are enormous. (2]||] There- 
fore, there is every reason to seek fundamental causes of 
these deviations, which are general for all HTSC. 

Several models have been suggested. Ovchinnikov 
and Kresin Q focused attention on magnetic impurities, 
which, as they assumed, cause pair breaking and effec- 
tively suppress superconductivity near T c . The tendency 
to magnetic ordering at lower temperatures results in a 
lower spin-flip scattering amplitude, thus enhancing su- 



perconductivity. The presence of magnetic impurities is a 
common feature of HTSC, since current carriers in most 
of them are due to doping, which generates magnetic de- 
fects at the same time. 

Spivak and Zhou || studied the role of Landau quanti- 
zation combined with a random potential. The quantiza- 
tion leads to a higher density of states on Landau levels, 
whereas the random potential brings to the Fermi level 
Landau sub-levels with opposite spins at points close to 
one another in space. In this case, the random potential 
must satisfy two opposite conditions: its variation over 
the coherence length £ should be larger than the Zeeman 
splitting, on the other hand, scattering by this potential 
should not wipe away peaks in the density of states. The 
HTSC structure favors both these conditions: fluctua- 
tions in the concentration of dopants, which are at the 
same time scattering centers, should occur even in high 
quality crystals, but these scatterers and current carriers 
are separated in space. 

It is possible that there are more fundamental causes 
of the peculiar shape of £? c2 (T) curves that can be put 
down to an exotic nature of superconductivity in HTSC. 
One example is the "bipolaron" or, in a more general 
approach, the "boson" model of superconductivity sug- 
gested by Alexandrov and Mott. Jl(| The model assumes 
that pairs (charged bosons, e.g., bipolarons) are pre- 
formed, and the superconducting transition consists in 
Bose-condensation of these pairs. In the presence of a 
random potential, the curve of £> c2 (T) has a positive cur- 
vature. The conventional superconductivity in a Fermi 
liquid can transform to the boson superconductivity if 
the electron-phonon coupling is strong and the carrier 
density is low. Again, HTSC materials are good can- 
didates for realization of such a scenario. Their carrier 
concentration is lower than in conventional metals and 
drops further with decreasing doping level, whereas the 
coupling constant A > 1. 
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Abrikosov suggested for HTSC a model whose cen- 
tral component is a saddle-like singularity in the electron 
spectrum. This model predicts, in particular, a posi- 
tive curvature of the B c2 (T) curve pj because the prob- 
lem becomes effectively one-dimensional due to the sad- 
dle point; as a result, the magnetic field's capability of 
destroying superconductivity is limited considerably. In 
the absence of the paramagnetic limit, the model yields 
the divergent function B c2 (T), but if the paramagnetic 
limit is taken into account, the critical field is limited to 
a finite value. 

The experimental data accumulated over recent years 
are insufficient for making an ultimate choice of one of 
these model. Further research is needed, and the present 
paper is a step in this direction. We present an investi- 
gation of the effect of a magnetic field on the resistivity 
of YBa2Cu3 06+a: single crystals at doping levels below 
the optimal one. The aim of this work was to measure 
the temperature dependence of B c2 in this material at x 
such that T c < 30 K and derive from these data changes 
in parameters that control B c2 when T — > 0. 

The paper is organized as follows. Section || presents 
basic theoretical concepts concerning the superconduc- 
tor phase diagram in a magnetic field and the behavior 
of conductivity around the superconducting transition 
point; they are essential in the analysis of experimen- 
tal data. Section III describes sample fabrication tech- 



niques and experimental procedures. Sec. IV reports on 
experimental results. The curves p(T) and their evo- 
lution induced by the magnetic field are discussed in 
Sec. IV A. The derivation of B c2 (T) from resistance- 



versus-temperature data for HTSC has remained a con- 
troversial issue, p^,p^[ t herefore this topic is given special 
treatment in Sec. [VE. Since the transition broaden- 



ing induced by magnetic field is insignificant, qualitative 
conclusions concerning the behavior of B c2 (T) are not af- 
fected by the specific routine employed in determination 
of the superconducting transition point. Nonetheless, in 
determining B c2 {T) quantitatively, we analyzed the fluc- 
tuation conductivity in t he no rmal state as a function 
of temperature. Section [V C discusses B c2 (T) derived 
from experimental data: the curvature of B c2 (T) curves 
proved to be positive throughout the available temper- 
ature range, including the close neighborhood of T c ; no 
signs of saturation in the low-temperature range have 
been detected; the experimental data are compared with 
existing models. 



II. BASIC THEORETICAL CONCEPTS. 

A. Phase diagram 

The phase diagram of a type-II superconductor in the 
B-T plane in the mean-field approximation contains a 
Meissner region, where magnetic field is fully ejected 



from a sample, a mized state region, where a lcttice 
of Abrikosov's flux lines exists and a normal metal re- 
gion. These regions are separated by lines of second or- 
der phase transitions: B C \(T) between the Meissner and 
mixed phases and B c2 (T) between the mixed state and 
normal metal. 

Beyond the mean-field approximation, thermal fluctu- 
ations of the order parameter slightly change the phase 
diagram configuration. Now, it contains a region of "vor- 
tex liquid," where fluctuations change largely the order 
parameter phase (which can be interpreted in terms of 
free motion of Abrikosov's flux lines), and a region of 
critical fluctuations close to B c2 (T), where the order pa- 
rameter amplitude fluctuates and its mean value changes 
rapidly with the temperature or magnetic field intensity. 
There are superconducting fluctuations above B c2 (T) 
also, but their amplitude is small and decreases away 
from the line of B c2 {T). The phase transition to the su- 
perconducting state with a long-range order established 
occurs on the boundary between the vortex liquid and 
vortex lattice [melting line B m (T)], whereas the curve 
of B c2 (T) determined in the mean-field approximation 
defines the line of a crossover from the normal metal, 
where the order parameter fluctuation amplitude is low, 
to the vortex liquid, where the magnitude of the order 
parameter is almost unity. [^4|-^6[ 

In conventional superconductors, the regions of criti- 
cal fluctuations and vortex liquid are quite narrow and 
essentially unobservable. The melting line B m {T) coin- 
cides with B c2 (T), therefore, the mean-field approxima- 
tion adequately describes the phase diagram. In HTSC 
the situation is different. Owning to the high critical tem- 
perature, small coherence length, and high anisotropy, 
fluctuations play a more important part, and the vortex 
liquid phase occupies a considerable region of the phase 
diagram, so B m and B c2 are separated. Since fluctuations 
broaden features of field dependencies of transport and 
thermodynamic properties at point B c2 , it is most diffi- 
cult to determine this point in experiment. Nonetheless, 
the value of B c2 (T) is still very important since this is the 
parameter that controls the behavior of thermodynamic 
quantities in the region far from the line of transition, 
where the mean-field approximation is valid. 

In materials with strong pinning, the phase diagram 
is further modified: the pinning destroys the order in 
the vortex lattice and transforms it to a vortex glass. 
The melting line is replaced by the "irreversibility line" 
B*(T), above which vortices are depinned by thermal 
fluctuations and move freely even at very low current den- 
sities, which results in a finite resistivity and reversible dc 
magnetization. Below B*(T) vortices are pinned in the 
low current limit, and the magnetization curve shows a 
hysteresis. 
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B. Resistive transition 

In high-temperature superconductors with optimal 
doping, curves of p(T)\s form a fan with a common tran- 
sition onset point, so the positions of the transition on- 
set are almost independent of the magnetic field. @Jl7| 
The drop in the resistivity around the transition onset 
is controlled by the contribution of superconductive fluc- 
tuations to the conductivity. The characteristic field of 
fluctuation suppression is B c2 , hence the shift of the tran- 
sition onset should follow the function B C 2(T). On the 
low-tempcrature side, the resistivity should vanish when 
the vortex motion is frozen. Qualitatively, the line on the 
B—T phase diagram where the vortex mobility becomes 
significant is the "irreversibility line" B*(T). Thus, the 
resistive transition is confined by the lines B c i{T) and 
B*(T) and is associated with the vortex liquid region on 
the phase diagram so that the fan-like appearance of re- 
sistivity curves is due to broadening of this region with 
the magnetic field while the line of B c i (T) is almost ver- 
tical. 

The breadth of the vortex liquid region, hence the tran- 
sition width, is determined by the relation between pin- 
ning and fluctuations. The vortex depinning is favored 
by the small coherence length £, high temperatures, and 
weak coupling between neighboring superconducting lay- 
ers of Cu02, i.e., by the high anisotropy. Variations in 
the doping level (carrier density n) to both sides from the 
optimal doping n op t, lead to lower T c and larger £. On 
the other hand, the anisotropy is stronger at lower doping 
and weaker at higher doping levels. The resistivity curves 
of overdoped HTSC samples with high carrier densities 
and low anisotropy are similar to those of conventional 
superconductors with strong pinning. P,p"8| 

The difference between over- and underdoped states 
was demonstrated by comparing L^-zSr^CuCU samples 
with different x. |l8| Whereas a magnetic field of B = 8 T 
broadened by 15-20 K the resistive transition in an un- 
derdoped sample with x = 0.08 and T c « 30 K, the tran- 
sition curve in an overdoped sample with x — 0.20 and 
approximately the same T c was shifted by magnetic field 
without changing its shape. []l8f This observation was 
confirmed by other researchers, |l9| , po|l who also reported 
that decreasing the oxygen content in YE^CuaOg+z thin 
films and single crystals considerably enhances effects 
originated from vortex motion, in particular, increases 
transition broadening in the magnetic field. All these ex- 
periments, however, used samples with T c > 40 K, and 
it remained unclear whether this tendency should persist 
in the range of low T c . 

There is an alternative interpretation of the resis- 
tive transition in cuprates, which attributes most of the 
change in the resistivity to a phase transition between the 
vortex liquid and vortex lattice (vortex glass) at B m (T). 
@,fn]|22) In this case, the resistive transition is decom- 



posed into a resistivity jump on the B m (T) line [well 
below B C 2(T)] and a crossover on line B C 2(T), |2l| which 
can produce only slight changes in resistivity. 

The high conductivity in the normal state of overdoped 
cuprates might in fact mask the transition from the nor- 
mal to vortex liquid state. Q But changes in transport 
characteristics around B C 2 are evident even in high qual- 
ity YBa2Cu306+x crystals with optimal doping and very 



weak pinning. |23| They should be the much more no- 
table in underdoped samples, whose conductivity in the 
normal state is essentially lower. 



III. EXPERIMENTAL 

YBa 2 Cu306+j; single crystals were grown by slow 
cooling the melt containing 10.0 to 11.4 wt.% of 
YBa 2 Cu306+a; and eutectic mixture of 0.28 BaO and 0.72 
CuO as a solvent with subsequent decanting of residual 
flux. For our experiments, we selected single crystals 
without visible signs of block structure and shaped as 
plates 20 to 40 thick with areas of several square mil- 
limeters. After oxygenating at 500°C, they had T c «90- 
92 K and fairly narrow resistive superconducting transi- 
tions with AT < 1 K. 

In YBa2Cu306+iE, current carriers (holes) are gener- 
ated in Cu02 planes as a result of capturing electrons 
in layers of CuO^ chains. The hole density depends on 
the oxygen content x and configuration of oxygen atoms 
in chains in CuO^ layers. Consequently, the carrier den- 
sity in YBa2Cu306+a; (along with the superconducting 
transition temperature) can be varied by two methods: 
changing the oxygen content and varying its ordering in 
CuO x layers. 

The technique for changing the oxygen content is the 
high-temperature annealing, and it allows one to produce 
the whole range of states from antiferromagnetic insula- 
tor to optimally doped superconductor. The annealing 
temperature at a given partial pressure of oxygen con- 
trols the oxygen content in a crystal and is a convenient 
technological parameter in processing superconducting 
samples. In order to reduce the oxygen content to 
x = 0.37-0.47, we annealed crystals in air at 700-800°C 
and then quenched them in liquid nitrogen to prevent 
exchange of oxygen with the atmosphere during cooling. 

The second technique allows us to vary the carrier den- 
sity over a relatively narrow interval by changing the av- 
erage length of oxygen chains at constant x. |25|,|26| In 
chains of finite lengths, there are q oxygen atoms per 
q + 1 copper atoms, hence, one has (q + l)/q electrons per 
oxygen atom. For this reason, oxygen atoms in shorter 
chains are less efficient in capturing electrons from Cu02 
planes. The average chain length can change owing to the 
high diffusion mobility of oxygen in CuO^ layers at the 
room temperature and above. Longer chains have lower 
energy, but they contribute less to the entropy, which 
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TABLE I. Samples 
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makes them less preferable at high temperatures. The 
balance between these two factors determines the aver- 
age chain length in equilibrium (hence, the number of 
holes) as a function of temperature. The relaxation time 
strongly depends on temperature, so rapid cooling freezes 
the oxygen configuration, thus fixing the carrier density. 
In real experiments, we heated crystals to 120-140°C and 
then quenched them in liquid nitrogen. This procedure 
notably reduced the number of holes in the sample, hence 
lowered T c . After that samples could be stored in liquid 
nitrogen for indefinitely long times without any changes 
whatsoever. If a sample was exposed to the room tem- 
perature, the carrier concentration increased gradually 
owing to oxygen coagulation in longer chains. This aging 
process could be monitored continuously by measuring 
the sample resistance at a constant temperature and in- 
terrupted at any moment by cooling the sample, thus we 
could obtain any intermediate value of T c . The aging of a 
sample at the room temperature for several days returns 
it to its initial equilibrium state. Since all restructuring 
processes in the oxygen subsystem proceed at relatively 
low temperatures, this method allows one to obtain a se- 
quence of sample states with minimal differences in con- 
figurations of defects and pinning centers. 

All in all, we have studied three crystals at several car- 
rier densities in each. The sample parameters arc listed in 
Table Q. The different crystals are numbered 1 to 3, their 
states with different oxygen contents are labeled a and 
b, and the quenching states are referred to as quenched, 
intermediate, and aged. The ratio between resistances at 
the room temperature and 50 K, when the free path is 
largely controlled by defect scattering, is a characteristic 
of crystal purity. This parameter of sample 2 is a factor 
of about three higher than in samples 1 and 3. Parameter 
B* r will be discussed in Sec. IV C . 



We measured the resistance in the ab plane using a 
four-terminal circuit. Since YBa2Cu3 06+a; crystals with 
low oxygen contents are highly anisotropic, it is very im- 
portant that the current be uniformly distributed over 
the sample thickness, so that only one component of the 
resistivity tensor is measured. To this end, the current 
contacts were fabricated over the entire surfaces of two 
opposite crystal faces. The contacts were made by a sil- 



ver paste and fixed by annealing before all thermal ma- 
nipulations designed to vary the hole concentration. The 
resistance was measured by the standard technique using 
a nanovolt-range lock- in amplifier at 23 Hz. The probe 
current was weak enough to ensure the linear regime and 
avoid overheating even at the lowest temperatures. The 
uncertainty in the geometrical factor restricted the ac- 
curacy of absolute measurements of conductivity to 10— 
20%, nonetheless, note that the geometrical factor of each 
sample was the same in all conducting states. 

Most of experiments were performed in a cryostat with 
a 3 He pumping system, which allowed us to vary the tem- 
perature between 0.3-300 K. j2?j At temperatures of 0.3- 
1.2 K the sample was immersed in liquid 3 Hc, at higher 
temperatures it was in the 3 He atmosphere at a pres- 
sure of several torr serving as a heat-exchange gas. The 
temperature was measured by a carbon resistance ther- 
mometer calibrated by a reference platinum thermome- 
ter, 4 He vapor pressure, and cerium-magnesium nitrate 
in appropriate temperature ranges. The magnetic field 
of up to 8.25 T was applied along the c-axis. 

Sample 3b in the aged state with low T c was tested in 
a dilution refrigerator at temperatures down to 30 mK 
and magnetic fields of up to 14 T. 



IV. RESULTS 

A. Temperature dependence of resistivity 

In our experiments on samples with T > 30-35 K 
(sample 2b quenched anh aged), we record fans of p(T)\s 
curves, similar to those reported by other authors. p9 j20 
In samples with lower T c , the effect of magnetic field 
on the resistive transition is radically different, and in 
this publication we concentrate on these effects, namely, 
the behavior of YBa2Cu3 06+a; samples in states with 
T c < 30 K (samples la, 2a, 3a, and 3b in all quench- 
ing states). In these samples, magnetic field shifts the 
transition without a notable broadening (Fig. ^), which 
indicates that the effect of vortex motion on the shape of 
transition curve no longer dominates. Nonetheless, the 
shape of the transition curve is affected by the magnetic 
field, and one can see on curves of temperature deriva- 
tives dp/dT plotted on the right of Fig. [l] that these 
changes are nonmonotonic. Since the normal state re- 
sistivity is almost constant with temperature, the peak 
amplitude on the derivative curve is inversely propor- 
tional to the resistive transition width. These graphs 
clearly show that, irrespective of T c (<30 K), the transi- 
tion width is maximum at about 13-14 K. If the zero-field 
T c is higher, the transition first shifts to lower tempera- 
tures with magnetic field and broadens (Fig. |l]a). Then, 
below 13-14 K, the transition narrows concurrently with 
its shift to lower temperatures. If T c is initially lower than 
13-14 K (Fig. [l]b and |l|c), the transition is narrowed by 
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magnetic field concurrently with its shift to lower tem- 
peratures from the start, and the slope of the transition 
curve in magnetic field becomes steeper than at zero field. 

The comparison between samples la and 2 a demon- 
strates that the nonmonotonic change in the transition 
width with magnetic field is a reproducible property and 
is little affected by the sample quality. The superconduct- 
ing transition temperatures of these two samples were 
driven to one value by annealing (Fig. |l|b and |l]c), but 
their parameters in the normal state were notably differ- 
ent. Sample 2a contained less impurities and structural 
defects, as a result, its resistivity around T c was twice as 
small (Fig. [l]), it dropped more rapidly in the process of 
cooling from the room temperature to 50 K (Table |) and 
showed a smaller increase in the range of lower temper- 
atures. Nonetheless, irrespective of all these differences, 
both the transition shift rate in magnetic field and the 
evolution of transition curves of these samples are sim- 
ilar. Narrowing of the resistive transition in an under- 
doped YBa2Cu3 06+a; with increasing magnetic field in 
this temperature range was detected by Scidlcr et. at. 
p5[ but, since their measurements were presented in a 
different form, it is difficult to compare them directly 
with our results. 

Such a behavior of transition curves is observed for all 
samples with T c > 6 K. In states with lower transition 
temperatures, we were not able to achieve sufficiently 
narrow transitions at zero magnetic field to measure T c 
and transition width. Therefore, the measurement data 
for sample 3b will be given and discussed separately in 



Sec. [VC. 



B. Derivation of B C 2(T) from resistance 
measurements 

The absence of the notable transition broadening in 
magnetic field in YBa2Cu3 06+a; samples with low T c in- 
dicates that, unlike samples with T c > 30-35 K, they 
have a narrower region of the "vortex liquid" on the phase 
diagram. The transition width, however, is not so small 
that it could be neglected in determining B C 2(T). Since 
the point B C 2 is not marked by a sharp feature on curves 
of p(T), there is a good reason to determine this point 
by fitting a theoretical curve describing the crossover be- 
tween the normal metal and vortex liquid to the experi- 
mental data. In developing this approach, let us consider 
the sample conductivity as a sum of the normal and fluc- 
tuation components: <j(T) — a n (T) + <jfi(T). 

The fluctuation conductivity <jfi in quasi- two- 
dimensional systems in zero field is usually described by 
the Lawrence-Doniach formula: 



13 K 




10 20 30 
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FIG. 1. Conductivity p a b (on the left) and its derivative (on 
the right) at various magnetic fields aligned with the c-axis. 

(a) Sample la in the aged state; applied fields (from right to 
left): 0, 0.06, 0.12, 0.23, 0.35, 0.6, 0.8, 1.2, 1.6, 2.2, 3.0, 3.8, 
4.6, 5.5, 6.7, and 8.2 T. 

(b) Sample la in the intermediate state; applied fields: 0, 0.06, 
0.12, 0.23, 0.35, 0.6, 0.8, 1.2, 1.6. 2.2, 3.0, 3.8. 4.6, 5.5, 6.7. 
and 8.2 T. 

(c) Sample 2a in the aged slate; applied fields: 0.12, 0.23, 0.5, 
0.8, 1.2, 1.6, 2.2, 3.0, 3.8, 5.5, 6.7, and 8.2 T. 
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where d is the interplane separation. Friedman et. 
al. pgj] show that, even in analyzing optimally doped 
YBa2Cu306+a;, crystals with the resistivity anisotropy no 
higher than 30-100, one can neglect the factor in brackets 
which takes into account effects of the third dimension 
and use Aslamazov-Larkin's expression for two dimen- 
sions: 

In oxygen deficient crystals, the anisotropy is up to (5- 
10) TO 3 , Ij therefore Eq.(§) is a fortiori valid through- 
out the temperature range in question, except the neigh- 
borhood of T c . 

There is no consistent theoretical description of 
o~fi(T, B) in nonzero magnetic field for arbitrary B C 2(T). 
Ullah and Dorsey analyzed a/i in a system with 
strong fluctuations in magnetic field and suggested a scal- 
ing expression for the fluctuation conductivity, which is 
often used in describing the resistive transition and deter- 
mining B C 2(T) of cuprate superconductors. p9|-j3l| Since 
their approach is based on the mean-field approximation 
and assumes a linear dependence B C 2(T) near T c , it does 
not apply when B c2 (T) is strongly nonlinear. (It will 
be shown below that this is the case in our samples.) 
Nonetheless, in the region well above T C (B) (es ^ 0.1), 
where Gaussian fluctuations dominate, a formula similar 
to that suggested by Aslamazov and Larkin can be used: 

cr/zcxe^ 1 , es=lii5^y (3) 

in both zero and finite magnetic fields (see Ref. jl6| and 
references therein). Here T C (B) is the functional inverse 
of B C 2(T). This formula also assumes, generally speak- 
ing, a linear dependence B C 2 (T) , but a possible change in 
the exponent of this function should lead only to a small 
systematic shift of the resulting curve T C (B). 

In contrast to the case of optimal doping, the nor- 
mal conductivity in our samples is low, of order of e 2 /Tid 
(Fig. |l|) , if d is assumed to be of order of the lattice con- 
stant, 11.7 A. Simple estimates based on the Aslamazov- 
Larkin formula (||) with a reasonable value of d indicate 
that the contribution of fluctuations, o~fi, should be sev- 
eral percent of o~ n even at es ^ 0.5. This makes de- 
termination of o~ n (T) more difficult. The difficulties are 
exacerbated by the fact that the normal state resistivity 
has a minimum in the region of 30-40 K and increases 
at lower temperatures. Therefore, we decide to select a 
priori the function a n {T) with several fitting parameters. 
The fitting to experimental data is performed by varying 
all parameters in both o/;(T) and a n {T). j2i| This pro- 
cedure could hardly produce sensible results if each curve 
p{T) were described by a different set of parameters. For- 
tunately, the magnetoresistance of YBa 2 Cu306+ I crys- 
tals in the discussed region of fields and temperatures is 
negligible in the normal state, i.e., the shape of o~ n (T) is 
constant with the magnetic field. 




FIG. 2. Characteristic points of the superconducting tran- 
sition in sample 2a plotted in (a) B-T and (b) p-T planes (the 
resistivity is normalized to function p n (T) in Eq.(^), which 
was used in determination of the fluctuation conductivity): 
(open circles) "irreversibility line", p = 5 /iSJcm; (squares) 
peak of derivative dp/dT; (full circles) _B C 2(T); (triangles) 
"transition onset" dajijdT = 10 2 Q~ 1 cm~ 1 /K. 

Our previous investigations of YBa2Cu30g+a: single 
crystals near the boundary of the superconducting region 
of the phase diagram |?2| revealed that the normal resis- 
tivity of such samples at T < 20 K is well described by 
a logarithmic function. In a broader temperature range 
(0.5 K< T < 150 K) the conductivity is very closely 
described by the empirical function 

a n {T) = p- 1 = [a~p\ogT + 1 T]- 1 (4) 

This function with three fitting parameters is used in 
processing our experimental data. 

By approximating the conductivity in zero magnetic 
field by a sum of a n (T) from Eq. (|4j) and er/;(T) from 
Eq.(||), T c and d being fitting parameters, along with 
a, (3 and 7, we obtain reasonable values d = 8-15 A, 
which are in fair agreement with the YBa2Cu306+a; lat- 
tice constant along the c-axis. This indicates that the 
Aslamazov-Larkin formula yields a correct estimate of 
the fluctuation conductivity in Cu02 layers and its ap- 
plication is justified. The normal conductivity is fitted 
so as to obtain the best approximation of the fluctua- 
tion conductivity throughout the range of magnetic field. 
Nonetheless, the uncertainty in the normal resistivity was 
quite considerable. It turned out, however, that calcula- 
tions of the transition temperature are little affected by 
admissible variations in er„(T). The resulting uncertain- 
ties in the transition temperature are shown in Fig. |^. 

This procedure enable us to derive B C 2(T) in the 
mean-field approximation from our measurements. Since 
the resulting curve of B C 2(T) is nonlinear and it casts 
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FIG. 3. Temperature dependence B C 2(T) in different 
states. Data for sample la in (empty triangles) quenched, 
(open circles) intermediate, and (squares) aged states; aged 
states of (crosses) sample 2a and (full inverted triangles) sam- 
ple 3a. 

doubt on the applicability of Eq.Q, we deem it neces- 
sary to demonstrate that, on the qualitative level, the 
shape of the B c2 (T) curve is not affected by subtleties 
of the data processing, owing to the absence of consid- 
erable transition broadening. Figure ||a shows the curve 
of B c2 (T) for sample 2a, along with its other characteris- 
tic fields, namely, the "irreversibility line" determined at 
p = 5 jtxOcm, positions of the peak of derivative dp/dT, 
and the line of "transition onset," which was defined as 
a point were daji/dT = 10 2 O^cm^/K. These lines 
are plotted in the B-T diagram in Fig. and Fig. |b 
shows positions of these points on the transition curves. 
(It is noteworthy that the values of B c2 are fairly close 
to those which would be obtained by defining the tran- 
sition point at a constant resistivity level p/pn = 0.8.) 
Figure |^a clearly shows that all curves in the B-T plane 
have positive curvature throughout the range of studied 
magnetic fields, including the region of low fields. This 
leads us to a conclusion that, even if the data processing 
procedure yields erroneous values of B C 2, the tempera- 
ture dependence of this parameter is qualitatively cor- 
rect. Our further analysis, however, will be based on the 
values derived from measurement data for the fluctuation 
conductivity. 

C. Universal temperature dependence of the upper 
critical field 

Measurements of B c2 (T) in three samples and five dif- 
ferent states (all the states of samples la and 2a and the 
aged state of sample 3a) are given in Fig. ||. It turned 




0.2 0.4 0.6 0.8 1.0 
T/T. 

FIG. 4. Dependencies B C 2 (T) for different samples reduced 
to the universal function b sc (t) using variables ([]). The nota- 
tion is the same as in Fig. ^. The inset plots the parameters 
Bsc (left-hand axis) and correlation length £o calculated by 
Eq. tj (right-hand axis). 

out that the curves for all the states can be brought to 
coincidence by varying the scales of the magnetic field 
and temperature, i.e., 

B c2 = B sc b sc (t), t = — , (5) 

J- c 

where B sc is the parameter characterizing the state 
and b sc (t) is a universal function (Fig. Function b sc (t) 
contains an arbitrary numerical factor. In Fig. ^ param- 
eter B sc is defined as B c2 at a specific reduced tempera- 
ture equal for all samples, namely, T c /2, i.e., the curves 
of B c2 (T) were brought to coincidence at two points, 
namely, at t = 1 and t = 0.5. The values of B sc for differ- 
ent states are listed in Table | and plotted in the inset to 
Fig. |J as a function of the zero-field transition tempera- 
ture. These points lie on one smooth curve, even though 
they are derived from measurements of the three differ- 
ent samples. The characteristic scale of magnetic field 
decreases (accordingly, the coherence length increases) 
with decreasing doping level more rapidly than T c , i.e., 
B sc is a super linear function of T c . This may be the main 
cause of the narrowing of the vortex-liquid region in the 
B-T phase diagram. As a result, YBa2Cu306+ K crystals 
with a high degree of underdoping with T c < 30 K do 
not display notable broadening of the resistive transition 
due to magnetic field. 

The function B c2 (T) was measured on sample 2b in a 
very narrow temperature range, T/T c > 0.9, owing to the 
limit on available magnetic fields. Its second derivative 
in this interval is also positive and all measurements of 
B c2 (T) can be fitted to function (0). But, since no data 
for lower temperatures are available and the expected 
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FIG. 5. Reduced resistivity of sample 3b in the aged state 
at temperatures of 50 and 36 mK as a function of magnetic 
field. The arrows indicate the difference between magnetic 
fields at which the curves cross the level p/pur = 0.8. 



critical fields are very high, the measurements of sample 
2b have not been analyzed in this context. 

Function b sc (t) is much different from bscsit)- First, it 
has no linear section near t = 1. This statement relies on 
Eq. (||), since for each B c2 (T) curve the limited precision 
allows one to draw a straight line of a small slope in 
the region within 1-2 K near T c , but if we consider the 
samples with higher T c , this linear region would be more 
narrow, and the slope of function b sc (t) at t = is smaller, 
which leads us to a conclusion that the universal curve 
has no linear section near t = 1. 

Second, b sc continues to rise as t — > 0. Figure |J shows 
this tendency in the region down to t = 0.1. In order 
to test the range of lower t, we investigated sample 3b 
with T ss 3 K at millikelvin temperatures. Its transi- 
tion curve is too wide to determine quantitatively T c and 
B C 2(T). Nonetheless, the measurements yield important 
qualitative information. Figure || shows the sample re- 
sistance versus magnetic field obtained at temperatures 
of 50 and 36 mK normalized to the resistance at a mag- 
netic field of 14 T. It is clear that a drop in temperature 
shifts the magnetoresistance curve to lower fields, i.e., 
B c2 (T) still grows with decreasing temperature even at 
T/T c ~ 0.01. We can obtain the following estimate: on 
the level pj p n = 0.8, which approximately corresponds to 
B c2 (T) according to Fig. |2|b, the magnetic field increases 
by 0.6 T; this yields a derivative of 40 T/K. Unfortu- 
nately, we cannot plot these points in Fig. Q for the lack 
of T c and B sc . 

This observation of B c2 (T) increasing even at very low 
temperatures is in accord with measurements of other 
materials, e.g., Tl2Ba2Cu06, where the upper critical 
field continues to grow at temperatures down to T/T c = 



0.001. 

Our data indicate that function B c2 (T) in under- 
doped YBa2Cu306+2; is nonlinear in the neighborhood 
of T c , and (dB/dT)\x c = 0. This conclusion contradicts 
most theoretical models based on the BCS model or the 
Ginzburg-Landau functional, which either predict a lin- 
ear behavior of this curve near T c or assume its existence 
a priori. This issue was not discussed in previous publi- 
cations of experimental investigations, but they all 
reported very low, if not zero, values of (dB/dT) at T c . 

The increase in the critical field owing to weakening 
of the spin-flip scattering predicted by Ovchinnikov and 
Kresin |^| should occur in the range of low temperatures, 
so it leaves the linearity of B c2 (T) near T c . essentially un- 
affected. The mechanism suggested by Spivak and Zhou 
H is effective only in high magnetic fields, where Landau 
quantization is significant, i.e., it also should not affect 
B c2 (T) near T c . Abrikoso derived B c2 (T) from the 
Ginzburg-Landau functional based on his model, which 
leads, naturally, to a linear dependence of B c2 in the first 
order in 1 — t. 

The nonlinearity of B c2 (T) near T c follows at present 
only from the model of bipolaron superconductivity 
Jl0| , p3]| which yields positive curvature of B c2 (T) for 
a charged Bose-liquid in a localizing potential, this 
throughout the entire temperature range. At temper- 
atures that are not overly low, the model predicts |33|] 

2 -| 3/2 

/ i' \ ' - i i \ • - 

B c2 (T) 



B*, ^ 



B* A 



2^1 



In) 



(6) 



Here £o is the correlation length and rii/2n character- 
izes the random potential. Equation (^) defines a uni- 
versal function in reduced variables without free param- 
eters. The only normalization parameter B* d corresponds 
to parameter B sc introduced in Eq. (|b|). It follows from 
Eq. (|^) that B^ — 0.68B SC . Comparison between our 
data and calculations by Eq. (|^) (Fig. ||a) shows excellent 
agreement in the region T/T c > 0.3. At lower reduced 
temperatures experimental points deviate from the the- 
oretical curve, but note that in this range we have only 
measurements of one state (sample 3a aged). 

The factor (1 — ni/2n) 1 / 2 in Eq. (^) is unknown, but, 
since neither in state 3a nor in state 3b have we detected 
a reentrant behavior of B c2 (T) predicted by Alexandrov, 
p3| it should be rather close to unity. Assuming this, we 
can derive from Eq. (||) the correlation length £o (Fig. 
right-hand axis in the inset). The length £o varies be- 
tween 70 and 300 A. The notable increase in the corre- 
lation length may be the main cause of the narrowing of 
the vortex liquid region on the B-T diagram. 

In the low-temperature region 0.1 < t < 0.6 the func- 
tion b sc (t) can be empirically described by the exponen- 
tial 
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FIG. 6. Function bsc(t) plotted in different coordinates: 
(a) the coordinates are selected in accordance with the bo- 
son model, Eq. (^); the inset shows the section close to 
t — T/T c = 1 on the extended scale; (h) semilogarilhmic 
coordinates; the dashed line follows Eq. (pi). 



In the low-temperature range T/T c < 0.3, experimen- 
tal points deviate from function (||) . On the other hand, 
measurements in the temperature interval between the 
lowest accessible values and t « 0.6 follow function (|?|). 
The combination of Eqs. (||) and (R) analytically de- 
scribes function b sc {t). 

However, the "universality" of function b sc (t) is lim- 
ited. We tested this function on our measurements 
of K .4Ba .6BiO3, Q and the experimental curve after 
renormalization of variables according to Eq. (^) was dif- 
ferent from function b sc (t) plotted in Fig. ||. 
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b sc = b cxp(-t/t ) (7) 

with parameters bo = 15 and to = 5.4 (Fig. |^b). Such 
an unusual temperature dependence in the region of low 
temperatures was detected previously |2^,|34| in measure- 
ments of the irreversibility line B* (T) . We suppose that 
the region of vortex liquid in phase diagrams of samples 
with low T c is a very narrow strip between B* and B C 2 1 
hence B* closely follows B C 2(T), especially at low tem- 
peratures. 

V. CONCLUSIONS 

Our investigation has supplemented the list of mate- 
rials displaying anomalous temperature dependence of 
the upper critical field B C 2{T) with underdoped cuprate 
YBa2Cu306+x- We have studied samples with differ- 
ent carrier concentrations and T c ranging between 6 and 
30 K. Throughout the studied temperature range, the 
curve of B C 2(T) for these samples have positive curva- 
ture and does not saturate at low temperatures. The 
curves for states with different T c can be brought to co- 
incidence in reduced coordinates T/T c and B / B SC (T C ). A 
fundamental feature of the universal function b sc (t) ob- 
tained in this manner is the tendency of its first deriva- 
tive dB C 2/dT to zero as T — > T c . Such a behavior can 
be interpreted at present only in terms of the model 
p0| , ^3[ treating the superconducting transition as Bose- 
condensation of preformed pairs. Other models designed 
to interpret the anomalous shape of the B C 2(T) curve 
predict a linear temperature dependence of B C 2 near T c . 
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